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Abstract 

Given a function / defined on a bounded polygonal domain Q C JR^ and a number A'^ > 0, we study 
the properties of the triangulation T/v that minimizes the distance between / and its interpolation 
J on the associated finite element space, over all triangulations of at most N elements. The error is 

studied in the W^''' semi-norm for f < p < oo, and we consider Lagrange finite elements of arbitrary 
polynomial order m—1. We establish sharp asymptotic error estimates as N ^ +oo when the optimal 
anisotropic triangulation is used. A similar problem has been studied in [4l llll l8l ll2U18| . but with the 
error measured in the norm. The extension of this analysis to the W^''' norm is required in order to 
' . ' match more closely the needs of numerical PDE analysis, and it is not straightforward. In particular, 

the meshes which satisfy the optimal error estimate are characterized by a metric describing the local 
aspect ratio of each triangle and by a geometric constraint on their maximal angle, a second feature 
• that does not appear for the error norm. Our analysis also provides with practical strategies 

'~~) for designing meshes such that the interpolation error satisfies the optimal estimate up to a fixed 

multiplicative constant. 

I ^ I Key words : anisotropic finite elements, W^''' norm, adaptive meshes, interpolation, nonlinear 

approximation. 
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> 

Introduction 

oo 

In finite element approximation, a usual distinction is between uniform and adaptive methods. In the 
latter, the elements defining the mesh may vary strongly in size and shape for a better adaptation to 
the local features of the approximated function /. Such procedures are used to improve the efficiency of 
] numerous numerical methods in scientific computing. This naturally raises the objective of characterizing 

. . and constructing an optimal mesh for a given function /. 

^ In this paper we consider a bounded bidimensional polygonal domain C M^, a fixed integer m > 2 

and an exponent 1 < p < cxd. For a given conforming triangulation T of 57 we denote by I™~^ the 

}_j standard interpolation operator on the space of Lagrange finite elements of degree to — 1 associated to 

T. A general objective is to study, for any / e C™(f7), the optimization problem 

^mf^^l|V(/-irV)llL.(n) (1) 

where the minimum is taken over all (possibly anisotropic) triangulations of cardinality < N. All the 
triangulations considered in this paper are assumed to be conforming: they have no hanging nodes, which 
implies that the interpolant I™~^ / is continuous and thus belongs to W^'°°{fl). The choice of the W^'P 
semi-norm appearing in the expression ([T]) is motivated by PDE analysis, e.g. elliptic equations in the 
case p — 2. Yet our paper in mainly a contribution to approximation theory, and aims at characterizing 
the approximation power of finite elements on anisotropic meshes. We obtain sharp estimates of the 
asymptotical behavior of the quantity ([T]) as — )■ oo, in terms of the TO-th derivatives of /. We also 
describe practical strategies for constructing meshes that behave similar to the optimal one, in the sense 
that they asymptotically satisfy as A^ — ^ c» this sharp error estimate up to a fixed multiplicative constant. 
Obtaining similar estimates and constructions in a non asymptotic setting remains an open question. 
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Estimates of a similar asymptotical nature were obtained in |4l [17] in the particular case of linear 
finite elements and with the error measured in the norm, instead of the VF^'P semi-norm. They have 
the form 

limsupfA^ min |1/ - < C ^\dei{d^f)\ (2) 

.1- . \ -LL I ^ AT 'I T n ( nA-\\ I 



iV-,.+oo ^ #{r)<N 



LP/(P+1)(0) 



which reveals that the convergence rate is governed by the quantity i/j det(d^/)|, which depends non- 
linearly on the Hessian d^/. This is heavily tied to the fact that we allow triangles with possibly highly 
anisotropic shape. The convergence estimate ([2]) has been extended to arbitrary approximation order in 
[15] , where the quantity governing the convergence rate for finite elements of arbitrary degree m — 1 was 
identified. This quantity depends nonlinearly on the m-th order derivative d™/. See also the book chapter 
|12j for an introduction to the subject of adaptive and anisotropic piecewise polynomial approximation. 

Main results and layout 

The Taylor development of a function / G C""(f2), close to a point z G 51, can be written under the form 

/(z + h) = y.,{K) + i^^ih) + o(|/;r ), (3) 

where h G is small, and where /ix and tt^ are polynomials which respectively belong to the spaces 

Pm_i := Span{x''?/' ; /c + ^ < to - 1} and ]H,„ Span{a;''?/' ; fc + ^ = to}. (4) 

For any triangle T, we denote by the local interpolation operator acting from C'^(r) onto Pm_i. 

For any continuous fonction v G C'^lT), the interpolating polynomial V^~^ v G Pm-i is defined by the 
conditions 

for all points 7 G T with barycentric coordinates in the set {0, • • • , 1}. If T is a sufficiently small 
triangle containing the point z, we thus have at least heuristically on T 

V(/-IJ'-V)^V(7r,-I™-V,), (5) 

since the Lagrange interpolation operator on the triangle T reproduces the elements of Pm_i. 

A key ingredient in this paper is the shape function Lm,p, which is defined by a shape optimization 
problem: for any tt G Hm, we define 

L„,p(7r) ^ mf^ ||V(^ - I^^i 7r)||i,(r), (6) 

where the infinium is taken over all triangles of area |r| = 1. The solution to this optimization problem 
thus describes the shape of the triangles of area 1 which are best adapted to the polynomial tt in the 
sense of minimizing the interpolation error measured in W^'^. If T is a triangulation of a domain J7, then 



I™^^ refers to the interpolation operator which coincides with I™^^ on each T G T. In view of (J5|, the 
optimization problem appearing in (|6| can be regarded as a "local" version of the "global" problem ([T]) 
of interest. 

The function Lm,p is the natural generalisation of the function Km p introduced in |18j for the study 
of optimal anisotropic triangulations in the sense of the interpolation error: for all vr G Hm 

Krn,p{'n) := I mf^ Ik - I™"-^ A\Lp{n)- 

Throughout this paper we denote by t G (0, 00) the exponent defined by 

1 TO — 1 1 , , 

T 2 p 

Consider two triangles T, T' and a polynomial tt G H^. If T is mapped onto T' by a transformation of 
the form z ^ az + (3, where a G 1R\ {0} and /3 G IR'^ (in other words the composition of a translation, an 
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homothety, and a central symmetry if a < 0), then recalling that tt is m-homogeneous one easily checks 
that 



ITI 



|V(^-I?,'-V)|| 



LP(T} 



|T'r^||V(7r-I™rV)|U.(T 



(T')- 



(8) 



Therefore ||V(7r — I™ ^7r)||Lp(T) ^ I?"! ^m.pl''') for any triangle T and any tt G Hm. Our asymptotic 
error estimate for the optimal triangulation is given by the following theorem. 

Theorem 1. For any bounded polygonal domain fl C IR^ , any function f € C"'(f2), and any 1 < p < oo, 

there exists a sequence of triangulations {Tn)n>No offl, with #(7Ar) < N, such that 



lim sup N 

N—foo 



1|v(/-i^-V)IU.(o) 



< 



L 



cTf 



(9) 



This theorem is the consequence of a sharper result, Theorem [2] which is given below. In the above 
estimate we slightly abuse notations by defining for each z 



L. 



7n,p 



d^fiz) 



:= L 



(10) 



where the polynomial tt^ G Hm is defined by ([s]). In other words we identify the collection d"'^f{z) of 
m-th derivatives of / at a given point z G to the corresponding term in the Taylor development of / 
close to z. Explicitly the right hand side of ([9| stands for 



L. 



m.p 



d'^f 



Lm,p{'^ z)''' dz 



The integer Nq appearing in Theorem [T] is independent of / and refers to the minimal cardinality of 
a conforming triangulation of Vt. An important feature of the estimate ([9| is the "lim sup". Recall that 
the upper limit of a sequence (un)n>No is defined by 

limsupuAr:= lim sup m„, 

and is in general stricly smaller than the supremum supjv>Ar(, un. It is still an open question to find an 

appropriate upper estimate of sup^>^^ W^if ^ ^Tn^ f)\\Lp{n) when optimally adapted anisotropic 

triangulations are used. 



We show in SjT] that a triangulation satisfies the optimal estimate of Theorem [T} up to a fixed multi- 
plicative constant, if it obeys the following four general principles: 

(i) The interpolation error should be evenly distributed on all triangles. 

(ii) The triangles should adopt locally a specific aspect ratio, dictated by the local value of d"^f. 

(iii) The largest angle of the triangles should be bounded away from 77 = 3.14159 . . . 

(iv) The triangulation T should be sufficiently refined in order to adapt to the local features of /. 
The third point (iii) is the main new ingredient of this paper compared to |18j , and is necessary for 

obtaining optimal W^'^ error estimates (but not for error estimates). Roughly speaking, two triangles 
having the same optimized aspect ratio imposed by (ii) may greatly differ in term of their largest angle, 
and the most acute triangle should be preferred when error is measured in W^'P rather than L^. The 
influence of large angles in mesh adaptation has already been studied in [6l [HI EH [9] . The heuristic 
guideline is that large angles should be avoided in general, since they lead to oscillations of the gradient 
of the interpolant. On the contrary, extremely thin triangles and very small angles can be necessary for 
optimal mesh adaptation. 

The shape function L„j,p plays an important role in our results, and we therefore devote f|2] to its 
study which is based on algebraic techniques. We obtain explicit minimizers, up to a fixed multiplicative 
constant, of the optimization problems which correspond to piecewise linear and piecewise quadratic finite 
element approximation. We also introduce, for arbitrary m > 2, explicit functions tt G i— Lm(7r) 
which are defined as the root of a polynomial in the coefficients of tt, and are uniformly equivalent to the 
shape function Lm,p, leading therefore to asymptotic error estimates similar to ([9]) up to multiplicative 
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constants. 



In order to illustrate the sharpness of Theorem [T] we introduce a slight restriction on sequences of 
triangulations, following an idea in a sequence {Tn)n>No of triangulations is said to be admissible if 
4I^{Tn) < N and s\xp j^;^ i^^{N ^ syxpj,^^^ diam(r)) < oo, in other words if 



sup diani(r) < CaN' 

TeTjv 



(11) 



for some constant Ca > independent of N. Here and below we denote the diameter of a set C IR 
by diam(i?) := sup{|x — y\; x,y £ E}. The following theorem shows that the estimate ([9| cannot be 
improved when we restrict our attention to admissible sequences of triangulations. It also shows that 
this class is reasonably large in the sense that ^ is ensured to hold up to small perturbation. 

Theorem 2. Let fl C IR^ be a bounded polygonal domain, let f G C"™(f2) and let 1 < p < oo. For any 
admissible sequence {Tn)n>No of triangulations offl, one has 



lim inf N 



-||v(/-i™-V)IU.(o) 



> 



(12) 



Furthermore, for all e > there exists an admissible sequence of triangulations {T§)n>No such that 



lim sup N 

JV-yoo 



|v(/-i^rV)ll 



Lp{n) 



< 



d'"f 



(13) 



Note that the sequences of triangulations (TJ^)Af>7Vo satisfy the admissibility condition (11) with a 
constant Cyi(e) which may grow to +cx) as e — > 0. Theorem [l] can be inferred from the estimate (|13|) 
proceeding as follows: for each N > Nq and each e > we define a real S{N, e) by the equality 



N- 



"l|v(/-i^rV)llL.(o)- 



m,p 



d"'f 



+ S{N,e). 



We next observe that for any fixed e > one has limsup^_j.oQ 6{N,£) < e. We may therefore construct, 
using a diagonal extraction procedure, a sequence {£n)n>No such that limsup^_j.o^ (S(iV, e^r) < 0. The 
sequence of triangulations {T^")n>No then clearly satisfies ([9]), which establishes Theorem [ij 



The proof of Theorem [2] is given in ^ The proof of the upper estimate ((13| involves the construction 
of an optimal mesh based on a patching strategy adapted from the one encountered in [1]. However, 
inspection of the proof reveals that this construction only becomes effective as the number of triangles 
N becomes very large. The construction described in ^|l] should therefore be preferred in practical appli- 
cations. 



Notations 

We denote by {u, v) the inner product of two vectors u,v & M^, and by |u| :— ^ {u, u) the euclidean norm 
of u. When g e LP{fl,TR^) is a vector valued function, we denote by ||.g||Lp(r2) the norm of a; n- \g{x)\ 
on ri, for instance in ([iJ and (|6|. 

We denote by M2 the set of all 2 x 2 real matrices, equipped with the spectral norm \\A\\ :— 
max|„|<i \ Au\. We denote by GL2 C M2 the group of invertible matrices, by SL2 C GL2 the special 
group of matrices of determinant 1, and by O2 C GL2 the group of orthogonal matrices. We denote 
by 5*2 C M2 the linear space of symmetric matrices, by 5® C 5*2 the subset of non- negative symmetric 
matrices, and by S*^ C 5® the subset of positive definite symmetric matrices. 

For any two symmetric matrices S, S' G 5*2 , we write S' < S" if and only if S" — S* S Sf . For any 
S € 5® (resp. S2) and any a > (resp. a G R) we denote by S" the symmetric matrix obtained by 
elevating the eigenvalues to the power a in a diagonalization of S. 

The greek letter tt always refers to an homogeneous polynomial tt € ]H,„, while the bold notation 77 
refers to the mathematical constant TT = 3.14159 . . . 
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1 Adaptive mesh generation 



This section describes some properties that are needed to ensure that a mesh T satisfies the optimal error 
estimate introduced in Theorem [l] up to a fixed multipHcative constant. We first introduce a description 
of the triangle based on some parameters adapted to our purposes. We then prescribe the behavior of 
these parameters in a triangulation tailored to the approximation of a given function /, and wc discuss 
the practical challenges encountered in the construction of such a mesh. 

1.1 Description of a triangle 

A triangle T C is determined by the collection of its three vertices ui, V2 and W3. In the context of 
adaptive mesh generation, we rather adopt the following parameters. The position of T is determined by 
its barycenter 

ZT := (wi + U2 + W3)/3. 



The area, the aspect ratio and the orientation of T are encoded in a symmetric positive definite matrix 
G S2 defined by the equality 

2 

l<i<3 

Last we shall introduce below a real S'(T) > 1 which is tied to the largest angle of T . 

If a triangle T' is mapped onto T by the change of coordinates z H> Az + zq, where A S GL2 and 
zo G M^, then one easily checks that 

Ut' = A^UtA. (14) 

From this point onwards we denote by Tcq the triangle of vertices (cos(2A;7r/3), sin(2fc7r/3))o<fc<2i which 
satisfies T-Lt^^ — Id. Combining this observation with (14 1, Proposition 5.1.3 in 20 establishes that for 
any triangle T 

|T|VdetHT = |Teq|, (15) 
and that there exists a rotation U ^ O2 (depending on T) such that the change of coordinates 

z ^ Un^z - zt) (16) 

maps T onto Tcq. Furthermore, as illustrated on Figure [l] (left), we have the inclusions 

{zt + u] u^Utu < 1/4} C T C £t := {zt + w; u^Utu < 1}. (17) 

The inclusion of two triangles T,T' therefore implies an inclusion of ellipses, hence an inequality on the 
associated symmetric matrices: 

r CT ^ AHt' > Ht- (18) 
We denote by p{T) g [1, cxd[ the measure of degeneracy of a triangle T, which is defined as follows: 

piT) ■■= ^J\\nT\\\\n^'\\. (19) 

Lemma 1.1. For each triangle T, one has 

, , diam(r)2 ^ ^ 

Proof. It follows from (15) that |T|/|req| = \\'H.t\V^\'H^^\\, and from ([l7| that IIH^'H = diam(£T)/2 < 
diam(T) < diam(i£'T). Combining this with (19) we obtain the announced result. □ 
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Figure 1: Two triangles and the associated ellipses associated by (17) (left). The interpolation points on 
the triangle Tq are aligned vertically (right). 



As illustrated on Figure [T] the fact that a triangle T is acute, or not, is not reflected on the ellipsoid 
£t or the matrix Ht- Since acute triangles play a priviledged role in finite element approximation we 
introduce the measure of sliverness S(T) of a triangle T, which is defined as follows 

S{T) := inf{||i/>|| Wtp^^W; ij) G GL2 s.t. the image of T by V is acute}. (20) 

(Where we refer to the image of T by the linear change of coordinates z ^ ijjz o\y M^.) The quantity 
S{T) can be regarded as the distance from T to the collection of acute simplices. It immediately follows 



from ([16|) and (jigf that 

1 < S{T) < p{T). (21) 

Then next expression gives an explicit expression of S{T) in terms of the largest angle of T, which shows 
that it is equivalent up to a multiplicative constant to the quantities previously introduced in [31 US] and 
referred to as (Tmin(T) and 1/ cos(0') respectively. 

Proposition 1.2. For any triangle T with largest interior angle 9, one has S(T) — max{l,tan |}. 

Proof. The result of this proposition is trivial if the triangle T is acute, we therefore assume that T is 
obtuse. We assume without loss of generality that the vertices of T are 0, au and f3v, where a,/3 > 0, 
u,v G R^, |u| = = 1 and {u,v) = cos9. Note that \u - v\ = 2sin(6'/2) and \u + v\ = 2cos(6'/2). 
Let ip G GL2 be such that the image of T by V' is acute. We thus have {ip{u),i/j{v)) > and therefore 
iV'(u) - V'(w)l < IV'(u) + It follows that 

ll^ll > X '^(^) + ^(")' > ^^^"(^/^) - tan ^ 

ll^^ll II''- II - + mu)^i;{v)\ - 2cos(0/2) ' ^^"^ 2' 

Therefore S{T) > tanf. On the other hand, let ip be defined by ip{u) = (0,1) and ip{v) = (1,0). 
Obviously the image of T by V' has one of its angles equal to 7r/2, and is therefore acute. One easily checks 
that ||?/)|| ll'i/'"^ II = tan(0/2) and therefore S{T) < tan |, which concludes the proof of this proposition. □ 

Most error estimates available in the literature, such as in [6i|T5], are designed to control the gradient 
interpolation error of a function on a triangle, by second or higher derivatives of the approximated 
function. Our purposes require a slight variant of these estimates, given in the next lemma, in which the 
gradient interpolation error is controlled by the gradient itself of the approximated function. 

Lemma 1.3. There exists a constant C = C{m) such that the following holds. For any triangle T and 
any f G W'^^'^iT), one has 

l|V(/ - I?^-l f)\\L^iT) < CS{T)\\^f\\L^^T)- 

Proof Let Tq be the triangle of vertices (0, 0), (1, 0) and (0, 1), and let g G W^'°°{Tq). We define g{x, y) := 
g{x, 0) and h{x, y) := g{x, y) —g{x, 0). Since g does not depend on y and since the Lagrange interpolation 
points on Tq are aligned vertically, as illustrated on Figure [ij the Lagrange interpolant I^~^ff does not 
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depend on y either. Futhermore, for all {x,y) € Tq, we have \h{x,y)\ — \ /J^^ ^{x,s)ds\ < \\^\\l°°(To)- 



Hence 



o xm — 1 
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Jo 


dy 


L~(To) 


dy 



< Coll l^-' hU^^To) < CoCi\\h\\Lo.(Ta) < C^Ci 



L°°(To) 



dy 



i°°(To) 



where the constants Co and Ci are the L°°(T()) norms of the operators g G Pm-i '"^ € I'm-27 and 
g eC°{To)^ I™-^ 3 e P™_i respectively. 

Let e be an edge vector of the triangle T. There exists an afhne change of coordinates ^' on H^, with 
linear part ip G GL2, such that 'l'(r) = Tq and ipe = eo, where eq = (0, 1) is the vertical edge vector of 
Tq. Noticing that 



(e, V o n = (e, V((I?|,-i 5) o *)) ^ (e^, (V g) o ^) 



dy 



o M/, 



we obtain 

rm— 1 



e,Virn5o^))l|Lo=(T) 



drm — 1 

dy 



<CoC\ 



L~(To) 



9?/ 



CoCi||(e,V(<?ovl/))||^^(^). (22) 



Applying this inequality to g = / o vji 1 we obtain that 

||(e, VI™-V)IIl=^(T) < CoCiIKe, V/)|Uo.(r), 



(23) 



for any edge vector e e {a, 6, c} of T. Defining the norm |u|t max{|(u, a)|/|a|, |(?;,6)|/|&|, Kw, c) |/|c|}, 
we obtain 

II |VI™-V|t ||L=e(T) < CoCi|||V/|t||l~(t)- (24) 

For any v G one has cos(6'/2)|'i;| < |ii|t < |w|, where 9 denotes the maximal angle of T. Indeed the 
upper inequality is trivial and the lower one is follows from the fact that at least one of the edge vectors 



makes an angle less than 6/2 with v. Combining this with (24), we obtain 

CoCi 



|vi™-V 



Since 9 > tt/S we have 
proof with C = 2CoCi + 1. 



< 2tan(6'/2) < 2S{T) according to Proposition 



1.2 



which concludes the 

□ 



Remark. The definition (20) of the measure of sliverness S{T), Propositions 1.2 and Lemma 1.3 have 
analogs for simplices of arbitrary dimension, see for instance the definition (3.87), Proposition 3.6.1 and 
Lemma 3.6.3 in 



1.2 Construction of a triangulation adapted to a function 

In order to identify the optimal shape of the triangles, our first step is to reformulate the optimization 
problem appearing in the definition ([6| of the shape function Lm,p. One easily checks using the invariance 
property ^ that for any tt € 



=inf{|rr^; Ts.t. iTp ? || V(7r - I^^^ ^)|Up(t) < 1}, 



(25) 



and the minimizers of this optimization problem and the original one are homothetic. 

The next lemma shows that the exponent p can be disregarded in the above expression, when one is 
only interested in minimizing it up to a fixed multiplicative constant. 

Lemma 1.4. There exists a constant c — c{m) > such that for any triangle T and any vr € Mm 

c\\v{iT-r^-^Tr)U^^T) < |rri;||v(^-i?^-V)|U.(^) < ||v(^-i™-V)|Uo.(^). 
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Proof. We denote by Tq an arbitrary but fixed triangle of area 1. For all .g€ L°°{To,'R^) we obtain using 
Jensen's inequality 

||.g||Li(To) < \\9\\lp{To) < I|5||l~(To)- (26) 

Furthermore, since all norms are equivalent on the finite dimensional space lP^_i, there exists a constant 
c = c(m) > such that c\\ fj,\\ (^Xa) — IImIIli(To) for all ^ S fm-i- Therefore 

c||mIIl~(To) < II/^IIli(To) < IImIIlp(To) < I1mI1l~(To), 

c||/^I1l~(t) < \T\^^\\^^\\L^{T) < \T\^'hA\LnT) < II/^IIl~(t), (27) 

where we used in the second line a change of variables from Tn to an arbitrary triangle T. We conclude 
the proof of this lemma by injecting = V(7r — I™^^ vr) in (27 1, where tt G IH,„ is arbitrary. □ 



For each fj, G W^^-i '^^ define 



l/.|l:=supM^, (28) 

z^o \zr 



and for each A G M2 we define /i o A G H„_]^ by o v4(z) := ii{Az). Note that o U\\ = for any 
[/ G O2, and that for A/ G 

||/ioM-5|| ^supMg^ where |z|m (29) 

hence for M, M' G 5^ 

M<M' ^ ||7roM"3|| > ||7roM'"5||. (30) 

The next lemma shows that, among the triangles having a prescribed aspect ratio, encoded by a matrix 
M, the smallest approximation error is achieved by the acute triangles (up to a fixed multiplicative 
constant) for which the measure of silverness S is minimal and equals one. 

Lemma 1.5. There exists a constant C — C{m) > 1 such that the following holds. For any triangle T 
and any tt G IH^n, one has denoting M :— "Ht 

C-i||(V7r)oM-5|| < ||V(7r-I^-V)||z,oo(T) < C5(r)||(V7r) o M-5||. (31) 
Proof. We may assume that the triangle T is centered at the origin of . We obtain combining Lemma 



1.3 and the inclusion T C St, see (17), 

|lV(7r - 7r)|U.o(i.) < C5(r)||V^|U=e(T) < ra(T)||V^|Uoo(£^) = C5(r)||(V7r) o A/-5|| 



which establishes the right part of (31). On the other hand, if C is sufficiently large, one has for all 

C < inf ll/i- i^||ioo(r ), 

since both quantities are norms on H^j_]^. Since zt = 0, the change of variables 2 M~iu'^z maps 
Teq onto T, where U G ©2 is defined by (16). Applying the above inequality to fj, o (M^^U^), where 
fj, G H^_]^ is arbitrary, we thus obtain 

C^^ll^o A/-5|| = C-^ll^o (M-5[/'^)|| < inf \\fio{M-^U^)-i^\\L^,T )= inf - j^||l~(t)- 



Choosing fi = Vtt and = V I™ tt we obtain the left part of (31 ), which concludes the proof. □ 
We introduce a variant Lm of the shape function, following an idea originally proposed in [8]: for each 

TT G Mrn 

Lmin) := inf{(detM)'^; M G 52^ s.t. ||(V7r) o M\\ < 1}. (32) 
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Lemma 1.6. There exists a constant C ~ C{m) > 1 such that for all vr G ]H„ 



Proof. Injecting (15), Lemma 1.4 and Lemma 1.5 into (25) we obtain 



inf{(|Teq|/|r|)^;Ts.t. \T\-^\\V{n -V^^' n)\\L.^T) < C} 

< inf{(detM)'^ ; M £ s.t. ||(V7r) o M\\ < 1} 

< inf{(|req|/|r|)"^;Ts.t. |Tr?||V(7r-IJ'-V)|Up(T) <c/C}, 

where we used the fact that for any M e S*^ there exists an acute triangle T such that Tir = M. Hence 

|rcq|'^L™,p(7r/C) < Lra{n) < \T,^\'-^L^.p{Cn/c), 

which concludes the proof since the shape function satisfies the homogeneity property Lm^p(A7r) = 
|A|L„jp(7r) for any A G M, as can be seen from the expression (|6|. □ 



Following a suggestion of the referees, we may compare Lemma [1.5| to other anisotropic error estimates 
that have been proposed in the literature. 



Remark 1.7. The specificities of the estimate (31) are the following: 



1. It only applies to an exactly polynomial function f = t: £ IHm, or f = tt + fi where fj, € Pm-i since 
the interpolation operator reproduces these elements: /i = I™^^ /i. This is a strong restriction, yet 
sufficient for our purposes, when compared for instance to IJj which applies to arbitrary functions 
in an adequate smoothness space. 

2. It is sharp ( up to the multiplicative constant C ) for any polynomial tt S Mm and any acute triangle 
T . The sharpness is lost in the case of a strongly obtuse triangle since the measure of sliverness 



S{T) only appears on the right of [31). A more elaborate estimate, which is restricted to the case 
m = 2 of linear interpolation but applies to general functions, is proposed in [T^ in the attempt 
(confirmed by some examples) to provide a sharp estimate for both acute and obtuse triangles. A 
detailed study of the interpolation error of a polynomial function, with respect to the shape of the 
triangle T , can also be found in |Pl I iOj / in the cases m = 2 or m = 3 respectively. 



3. The orientation and the scales of the triangle T are encoded in the matrix M = Hj 



e St- 



This 



is one of the major strengths of this estimate, because it leads to the optimization problem (32) 



posed on the set S*^ of symmetric positive definite matrices which can be addressed mathematically 
This problem is studied in using algebraic techniques which yield an explicit equivalent of the 
shape functions Lm,p for all m > 2, and an explicit "near minimizer" tt G Mm i— )■ AimiT^) € *S'^, 
for m G {2,3}, of the minimization problem (32) defining LmiT^)- This problem is also studied 
in Chapter 6 of 120^ using analytical techniques, which yield well behaved, although implicit, "near 
7W„(7r) of for all m > 2. 



mmimizers tt 



We now focus our attention on the global approximation of a function / G C""(ri) in the W^'P 
semi-norm. For that purpose we assume that a map z £ ^1 i-^ M{z) & S2 has been obtained which 
satisfies 



{detM{z)y-^^dz<Cl Lm{TTzVdz, and \\tTz o Miz)'^\\ < 1 for all z e 



(33) 



where the polynomial tt^ G H„i is defined by (|3|, the exponent r by ([t]), and Cl is a constant not too 
large. In other words we assume that the matrix M{z) is a minimizer, in an average sense and up to 
the sub-optimality constant Cl, of the optimization problem appearing in the definition (32 1 of Lm{T^z)- 
Such a map can be obtained by setting M{z) :— M.m{T^z) + <5Id, where Mm is described in Point 3 of 



Remark 1.7 and where (5 > is a positive constant introduced to avoid degeneracy problems. 
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Figure 2: Anisotropic mesh generation with a guaranteed upper bound on the maximal angle of the 
elements is generally tractable if one only requires some anisotropy close to the boundary of the domain, 
and tangentially to it (left). This is not any more the case if some anisotropy is required in the interior 
of the domain, at least with current software such as \22^ (right). 



We introduce a Riemannian metric H : ft S2 

H{z) := h-^ (det Miz)) (— Miz), 



(34) 



where /i > is a parameter. Mesh generation software such as [22] j see Figure [2] (right), and [23^ in three 
dimensions (see [17] for a more extensive list), are designed to produce a mesh T of £7 such that 



C'o^Hiz) <Ht< ClH{z) 



(35) 



for all r e T, 2 G T, where Cq > 1 is a constant not too large which reflects the quality of the adaptation 
of the mesh T to the metric H. In the expression (34) of the Riemannian metric H, the matrix M{z) 
is used to prescribe optimal aspects ratios (requirement ii in the introduction) for the triangles T £ T, 
and the scalar factor to equidistribute the interpolation errors among the elements of T (requirement 
i). Theoretical guarantees for such algorithms were established in [TBI [71 [2^ when the metric sufficiently 
regular. Unfortunately these results do not guarantee any property of the measure of silverness S{T) of 
the generated triangles T G T (requirement iii) . This generally forbids to achieve the optimal convergence 
estimate stated in Theorem [T] even up to a fixed multiplicative constant. The adaptation of the method 
presented here to this (suboptimal) context is described [T^ . 

Limited results exist nevertheless on anisotropic mesh generation with some control on the measure 
of silverness. Such a construction is (usually) possible if one only requires some anisotropy close to the 
boundary of the domain, and tangentially to it, see Figure [2] (left), which is the adequate behavior for 
the discretization of number of problems (such as the Poisson equation, or singularly perturbed reaction 
diffusion problems) as discussed in [2]. In the general case where some anisotropy is required in the 
interior of the domain, we refer to Theorem 6.1.2 in Chapter 6 of [20], on "quasi-acute triangulations" , 
which implies the following. Assume that is the periodic domain (M/K)^; the key is the absence of 
a boundary. Assume that the metric H satisfies the strong regularity properties required for this result, 
which are expressed under the form of Lipschitz regularity with respect to some distances on S2 and , 
and are satisfied if the parameter /i > in (34| is sufficiently small according to Lemma 6.5.12 in |20) . 
Then there exists a triangulation T which satisfies (351, for an absolute constant Cq independent of H, 



and in addition the following property: there exists a Co-refinement T' of the triangulation T (in other 
words each element of T contains at most Cg elements of T', and each element of T' is contained in an 
element of T) such that S{T') < Cq for all T' G T' (requirement iii in the introduction). 

Let {T,T') be a pair of triangulations satisfying the above conditions. If T{z) E T denotes the 
element containing a point z G f2, then we obtain using (15) and (35) 



#(r') < Co#(r) = Co 



dz 



< 



cql Jn 



[ Vd^^^^h-' [ 



(detX)- 



(36) 



Let z€^l,T'(zT' and T e T be such that zeT <zT. Using (18) and (35) we obtain H'j, > ^-^Ut > 
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{2Coy^H{z), hence 

||VK-I™-V3)|U.(j.,) < \T'\i\\V{7T,-l^r'7T,)h^^T') 

< C|req|?(detHT')^^^(T')ll(V7r,)oHj;,^|| 

< Ci{detH{z))-^\\{WTT,)oH{z)-^^\\ 
= Cih^\{WTr,)oM{z)-^\\ 

< Ci0, 



where Ci C\Tcq\p2rCQ . We used successively Jensen's inequahty in the first line, (15l and Lemma 



1.5 in the second line, S{T') < Co the definition ([?]) of r and (30) in the third line, (34) in the fourth line 



and ( 33 ) in the last line. In order to control the function /, instead of the polynomial tt^ which corresponds 
to its Taylor expansion, we may use Point i of Lemma |3.1[ proven in Sj3]below, which immediately implies 
in this context that 

ll|v(/ - rr'r' /)IIlp(t') - l|v(7r, - r^r' ^.)||lp(t')I < ^h^, 

where e{h) — as ft- — 0, and the function e : ]R\, — 1R+ only depends on /, m and A4. If ft, is sufficiently 
small (requirement iv in the introduction), then e{h) < 1 and therefore || V(/— 1™,~^ /)||lp(T') < (Ci+l)ft? 
for each T' e T'. It follows that 

#(r')"^l|v(/-i^rV)llL.(o) < #(r')^ max ||v(/-i™-V)IU.(T') 

< C2||L,„(d"//™!)llLno), 



where C2 = (Ci + 1){Cq /\Tcq\)^ C^. We used the definition ([7|) of the exponent r in the first line, (36) 
in the second line and ( |33[ ) in the last. As announced the triangulation T' satisfies the optimal error 
estimate of Theorem [T] up to the multiplicative constant C2. 

2 Study of the shape function 

This section is devoted to the close study of the shape function Lm,p, using algebraic techniques. Our 



approach is based on the variant Lm introduced in ( 32 ) , defined for all tt e Hm by 



L^{tt) := inf{(detM)'^; Me 5*+ s.t. || (Vtt) o Af^s || < 1}^ (37) 



and which is uniformly equivalent to Lm,p according to Lemma 1.6 



2.1 Explicit minimizers 



We describe the solution to the optimization problem appearing in (37) when m S {2,3}. The case 
m = 2, which corresponds to piecewise linear finite elements is already known, and discussed in detail in 
[TBI [3] and [2] for functions of more than two variables. In contrast the results in the case m = 3 are 
entirely new, although this case, which corresponds to piecewise quadratic elements, has already been 
discussed in [TU| . 

In order to present our results, we introduce some notation. For any homogeneous quadratic polyno- 
mial vr e IH2, TT = ax^ + 2bxy + cy^ , we define the symmetric matrix 



a b 
b c 



(38) 



For all vr e H2, tt — ax^ + 2bxy + cy^ , we define 



a b \^ , / + &^ ab + bc 
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For all vr € IH3, tt ~ ax^ + ibx^y + icxy^ + c?y'^, we define 

2 




(40) 

We say that a polynomial vr e is univariate if there exists A, a, /3 g II such that tt = A(Q:a; + jSyY^ . 

Proposition 2.1. i. // vr G _/f/2 *s nof univariate, then the matrix M 2 (t^) is the unique minimizer of 
the optimization problem appearing in (32). 



ii. The map tt G — !■ A^3(7r) is a near-minimizer 0/ t/ie problem (j'^ m f/ie following sense. If 
TT G if/3 is not univariate, then M.^{tt) is non- degenerate and ||(V7r) o A^3(7r)^2|| < (hence 
IKVtt) o (y27W3(7r))-5|| <l). Furth ermore there exists a constant C, independent o/tt, such that 



yjAetMzi^) < CLsin). (41) 



Proof According to ([29|), ||(V7r) o M-i\\ < K is equivalent to |V7r(z)| < K\z\'J(f^ for all z G M^, where 
TT G Hm, M G S2 and if > are arbitrary. 

We first consider an homogeneous polynomial tt G H2, which is not univariate. For all z G one 
has V7r(z) = 2[7r]z, and therefore |V7r(z)p = z"^ Ai2{'^)z. On the other hand consider M G S2 such 
that |V7r(z)p = z'^A42{'^)z < z^'^Mz. We thus have Ai2{T^) < M in the sense of symmetric matrices, 
which implies that det A^2(7r) < detM with equality if and only if A^2(7r) = M, since A42{tt) is positive 
definite. This concludes the proof of the first point. 

We now consider tt G IH3, which is again not univariate. In the sense of symmetric matrices, we have 



where we used the fact that the square root : S"® — > Sf is increasing. It follows that 

|V^(z)P = |9.7r(z)|2 + |9,7r(z)p < 2{z^ Msi7r)zf , 
hence M3{tt) satisfies the constraint ||(V7r) o A^3(7r)~5|| < \/2. Note that 

det7W3(7r) = 9y/ {a^ + 2b^ + c^){b^ + 2c^ + d^) - {ab + 2bc + cd)^ . (42) 



We postpone the proof of (41 1 to S2.2 right after (48), as we develop a general method for obtaining 
simple equivalents of the functions Lm- O 

Remark 2.2. It was proposed in |3 | to generate an anisotropic mesh T of a domain via the 
transport of a uniform mesh T' of an auxiliary domain fl' by a diffeomorphism F : fl ^ fl' . Without 
entering the details of this approach, we may describe one of its successes. Assume that one wishes to 
approximate a strongly convex function f G C'^(Jl), using linear finite elements. Define F(z) := 2V/(z) 
and fl' := F{fl) C . Consider a uniform mesh T' of il', of mesh size h > 0, and denote by T 
the collection of triangles obtained as follows: for each triangle T' G T', of vertices Vi,V2,V3, the set 
T contains the triangle T of vertices F^^ [vi] ^ F~^ (V2) , F~^ [v^) . If the parameter h > is sufficiently 
small, then T is a triangulation si nce F : £7 — > fi' is a diffeomorphism. Furthermore we obtain using (I4) 
that T is adapted in the sense of (35) to the metric iJ : — > S*^ defined by 

H{z) := h-^dFizfdFiz) ^ h-^A{d^f{z)f = h-^M2{T^z), (43) 



hence the optimal metric, see (34), for the approximation f in the W^'°° semi-norm! (The control of the 
measure of sliverness of the elements of T remains, however, an open problem with this method.) 

Let us finally mention that, although they are derived from the coefficients of tt, the maps tt i— > jM„i(7r) 
for m G {2, 3} are invariant under rotation, and therefore not tied to the chosen system of coordinate 
{x,y), as expressed by the following result. 
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Proposition 2.3. For any m G {2, 3}, any tt g Mm and any unitary matrix U G O2, one has 

Mra{^oU)^U^Mra{T^)U. 

Proof. Let TT e H2 and let [/ e Then, as announced, 

M2{t^ o [/) = 4[7r o Uf = 4(C/'^[7r][/)2 = [/T7W2(7r)t/. 
Let TT G E3 and denote by (uij)i<ij<2 the entries of the unitary matrix U . Then 

[d^{-KoU)] = uiiU'^[d^T:]U + ui2U^[dy-K]U and [dy{-KoU)] = U2iU'^[d.^'K]U + U22U'^[dyiT]U 
Hence, since U is unitary, 

[d^{TT0U)f + [dy{TT0U)]'' 

= {ul^+ul^)U^[d^nfU +{uiiUi2+U2lU22)U^m.T^][dy^] + [dyTT][d^n])U + {ul^ + ul2)U^ [OyTlfU 

Therefore 



which concludes the proof. □ 
2.2 Polynomial equivalents 

We introduce equivalents of the shape function tt L,„(7r) on Hm, which can be written in analytic 



form in terms of the coefficients of tt e ]H„j. As a starter we infer from Point i of Proposition 2.1 that 
for any tt £ H2 

L2(7r) = (detX2(7r))i = 2^1 det^]. 

For each r > 2 we denote by the space of homogeneous bivariate polynomials of degree r, as in 
Q for r — m, equipped with the norm 

II II 

\\m — sup——. 

For 1.1 E Mr and A e M2, we define the homogeneous polynomial n o A E by /i o A{z) := fi{Az), 
z £ M^. Observe that for M e 5^ one has 

ll^oAf-^IHsup"^^")" 



z^o \z\ 



M 



We now introduce a variant Kr '■= H,. — )■ ]R+ of the shape function L„j, which was first defined in [8], 
and later studied in |18j , in the context of optimal mesh adaptation for finite element approximation in 
the LP norm. (More precisely, due to different conventions, the function is tied to the function 
defined in [TH] by the equality Kr = TT^iKr ■) For each ^ G we define 

Krin) inf{(detM)5; Me s.t. ||/ioM-3|| < 1}. 

Observe that |V7rp :— {OxTt)'^ + {dyiiY E H2m-2 for each tt G H,„, and that clearly 

Lmin) = v/^2™-2(|V7r|2). (44) 

Given a pair of non negative functions Q and R on Hm, we write Q ^ i? if and only if there exists 
a constant C > 1 such that C^^Q < R < CQ holds uniformly on Mm- We sometimes slightly abuse 
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notations and write Qijr) ^ i?(7r). We say that a function Q is a polynomial on ]H,„ if there exists a 
polynomial P of m + 1 real variables such that for all oq , • • • , am G IR, 

We define degQ := degP, and we say that Q is homogeneous if P is homogeneous. 
The following equivalences were established in flF : for tt g H2 

K2{7t) ^ VldctMl, (45) 

and for vr G H3 

K^in) ^ ^|disc(7r)|, (46) 

where disc(7r) denotes the discriminant of the cubic polynomial tt, which is defined by 

disc(ax^ + Sbx^y + icxy^ + dy^) = 4(a6 - c^){hc - d?) - {ad - hcf . 

More generally for each r > 2 an explicit homogenous polynomial Q on is introduced in 18 , which 
satisfies 

Kr ^ yiQI , with d := deg Q. 



Combining this result with ( 44 ) we obtain an explicit equivalent of the functions L„j . 

Proposition 2.4. Let m > 2 and let Q be an homogeneous polynomial on if/2m-2 such that K2m-2 ~ 
•^/jOj, where d ~ degQ. Let be the polynomial on Mm defined by 

Q.iTT) :=Q(|V^n, 

then Lm ^ ^\fQl on Mm ■ 

This construction uses an equivalent \/[Qf of K2m-2 to produce an equivalent of L„i- Unfortunately, 
as m increases, the practical construction of Q becomes more involved and the degree d quickly raises. 
In the following theorem, we build an equivalent to L„j from an equivalent of Km-i instead of K2m-2i 
which is therefore simpler. 

Theorem 2.5. Let to > 3 and let Q be an homogeneous polynomial on Mm-i such that Km~i ^ y/lQl? 
where d — degQ. Let {Qk)o<k<r be the homogeneous polynomials of degree d on M„i-i x M„i-i such 
that for all u,v Cz IR and all 7^1,7:2 G Mm one has 

Q{u7r,+vn2)^ ^ uS'^-'^-gfe(^i, ^2), (47) 

a<k<d ^ ^ 

where (J^) :— jTfjTfey!- Q* polynomial defined for all tt G Mm by 

0<k<d ^ 

then Lm ^ ^\fQl on Mm ■ 

Proof. See Appendix [B] □ 



Using this construction and (45) we obtain an equivalent of L3 as follows. Let tti = ax^ + 2bxy + cy^ 
and TT2 ~ a'x^ + 2b' xy + c'y^ be two elements of H2. We obtain 

det([M7ri + v'K2]) — {ua + vd){uc + vc) — [ub + vb'Y 

= u^{ac - b^) + uv{ac' + a'c - 266') + v'^{a'c' - b'^). 
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Applying the construction of Theorem 2.5 to tt = ax^ + 3bx^y + icxy^ + dy^ € H3 we obtain 

L3(7r) - 3^(ac - b^Y + {ad - hcf /2 + {bd - d^Y . (48) 

Remarking that 

2[(ac - h'f + {ad - bcf /2 + {bd - Yf] = {a^ + 29 + c^){h^ + 2c^ + d^) - {ab + 26c + cdf , 



and using equation (42) we obtain that i3(7r) ^ det J\A^{-k) . This point concludes the proof of Propo- 



sition 2.1 and thus establishes that the map A^a defined in (40) can be used for optimal mesh adaptation 



for quadratic finite elements. 



Using (46) we similarly obtain an equivalent of L4. Denoting tt — ax + Abx y + Qcx y + A^dxy + ey . 

Li{TTf - {ib^c^ - 4ac3 - Ab^d + Qabcd - a^d^f 

+ (2&c3 - QaYd + iabd^ - 4b^e + Gabce - 2a^def/A 

+ (3c'' - Qbc^d + 86^^^ - Qacd^ - Gb^ce + 6ac^e + 2abde - a^e^f/6 

+ {2c^d - 4ad^ - 6bYe + Ab'^de + Qacde - 2abe^fl'i 

+ (Bc^d^ _ 46^3 _ 4^3g ^ g^^^g _ b^e^f. 

The following proposition identifies the polynomials tt G Hm for which L,„(7r) = 0, and therefore the 
values of d™/ for which anisotropic mesh adaptation may lead to super- convergence. 

Proposition 2.6. Let m> 2 and let t.^ ■— [^^^J ■ Then for all tt G JH^n, 

Lm{'^) = i/ and only if n — {ax + fSyY"^?: for some a,/3 Ft and tt G lH„i^tm- (49) 

Proof. It was established in [TB] that if2m-2(2'*) = if and only if tt* G H2,n-2 has a linear factor of 
multiplicity m. We therefore obtain, using (44), that im(7r) = if and only |V7rp is a multiple of l™, 



where I is of the form I — ax + /3y. 

Let us first assume that |V7rp = {dxT^Y + {dyT^Y has such a form. Clearly (9^''')^ and {dyirY are 
both multiples of l"\ Therefore d^ir and 9j,7r are multiples of where s is an integer such that 2s > m, 
hence s > — 1. We therefore have 

dxT^ — l'''^! and dyir = Z*7r2 where 1^1,1:2 G H,„_i_s. 

Recalling that I — ax + Py we obtain 

= dlxTT - dly-K = l^dyTTl - 8x712) + sF-\l3Tri - aTT2), 

hence /Jtti — aTr2 is a multiple of I. Since tt is homogenous of degree m it obeys the Euler identity 
TO7r(z) — {z, Vtt{z)) for all z = {x, y) G H^. Assuming without loss of generality that a 7^ 0, we therefore 
obtain 

m7r(a;, y) = P(x7ri -I- yTT2) = \ {ax + (3y)— + —{an2 - /Stti) ) 

V a a ) 

which shows that tt is a multiple of , hence of /*". 

Conversely if tt is a multiple of then d^T^ and dyi: are both multiples of Since 2{tm — 1) > m 

the polynomial [Vvrp is a multiple of which concludes the proof. o 



3 Proof of the main result 

This section is devoted to the proof of Theorem [2j We thus consider a fixed bounded polygonal domain 
C H^, an integer m > 2, an exponent 1 < p < oo and a function / G C™(f2). 
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The Taylor development of / close to a point z € Q is given by two polynomials fiz G IPm-i and 
Ti'z G K^: f{z + h) = Hzih) + T^zih) + o(|/i|™), where h € is small. The Taylor development of the 
function V/ : fl — > close to a point z e is obtained by derivation of the previous one: 

V/(z + h)^ Wfizih) + Viizih) + od/il™-!), 

and the corresponding Taylor formula reads as follows: for any z, h such that [z, z + /i] C 

Vf{z + h)^Wfiz{h) + {m~l) f W7Tz+th{h){l-tr-^dt. (50) 

Jt=o 

We denote by oj the modulus of continuity of the function z £ Q i-^ VtTx e E^j_]^: for each r > 

w(r) := sup {II Vtt^ - Vtt^' II; z,z' e f7 s.t. |z-z'|<r}, (51) 



where we recah that \\n\\ = sup{|Ai(z)|/|z|'"-^ z 0} for all fj. G H^_i, see (28) 



ation error of /. 



Our first lemma is an estimation on a single triangle of the gradient interpo . 
Lemma 3.1. There exists a constant C* = Ci,{f,m) such that the following holds, 
i. For any triangle T d fl, and any z ^ T one has 

|||V(/ - 1"t'-' /)||lp(t) - ||V(^, - 7rz)\\L.(T)\ < aw(diam(T))|r|^p(r)"^ (52) 
ii. For any triangle T d one has 

l|V(/ - r^-' /)||lp(t) < a|T|^ diam(T)"-i5(r). (53) 

Proof. The point z is fixed throughout this proof. We define a function g G C™(r2) by the equality 

g{z + h) := f{z + h) - ^izih) - TTz{h), 
for all h such that z + ft, G fi. Wc have 

|||V(/ - ir ' f)\\LHT) - l|V(vr, - lTz)\\LnT)\ < I1V((/ - TTz) - l^-\f - TTzML^iT) 

= l|V(5 - I?"^5)IIlp(T) 

< |r|i||v(.9-i™-i5)IU==(T) 

< Cp{T)\T\-p\\Wg\\L^(T), 

where we used the reverse triangle inequality in the first line, the translation invariance® and the 



1.3 



in the last 



equality = I™ ^ in the second line, Jensen's inequality in the third line, and Lemma 
line combined with the inequality 1 < S{T) < p{T), see (21 1. We have 

V5(z + h) - V/(z + h)- \/pz{h) - \7TTz{h) = (m - 1) / {\/TTz+th{h) - V^,(ft))(l - ty^'-^dt, 

Jt=o 

therefore ||V5(z + h)\\ < \h\"'-'^oj{h). Furthermore if z G T and z + /i G T, then \h\'^ < diam(T)2 < 
(4/|Teq|)|T|p(T) according to Lemma [l.l[ This concludes the proof of Point i, provided that C* > 

C(4/|Teq|)^. 

We now turn to the proof of Point ii, and for that purpose we consider a fixed point z G T. Changing 
our previous notation we denote by g G C™{il) the function defined by g{z + h) :— f{z + h) — lJ.z{h) for 
all h such that z + /i G We obtain 

l|V(/ - V^~^ /)|U.(T) = ||V(g - l^-' g)\\L.(T) < |T|^ ||V(.g - r^-' .9)||l==(t) < C\T\Hv g^^^r). 

where we used successively that the interpolation operator reproduces the elements of Pm_i, Jensen's 
inequality, and Lemma 1.3 On the other hand we obtain using (50). 

V.g(z + h) - V/(z + /i) - V/i,(ft) = (m - 1) ( V^z+Mil - tf-^dt, 

Jt=o 

hence |Vg(z + h)\ < |ft|™"i supfUVTr^H; z e fl}. If z G T and z + ft G T, then |ft| < diam(r), which 
concludes the proof of (ISSl) provided that C+ > C sup{|| Vtt^H ; z G fi}. □ 
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3.1 The lower error estimate (12) 



Under the hypotheses of Lemma 3.1 and recaUing from (|8|) that |lV(7r — V^^ t^)\\lp{t) > l^^l ^ Lm.p{Tr), 
we obtain 

(54) 



I|V(/-IJ'-V)IIlp(t) > \T\^ (L™,p(^,)-ac^(diamrMT)^) 



We consider an admissible sequence of triangulations {Tn)n>No- For all N > Nq, all T ^ Tn and all 
z e T, we define (j)N{z) ■— \T\ and 



iPn{z) (i™^p(7r,) - ac^(diam(T))p(T)^ 

and the exponents Pi — ^ and P2 ~ (m-Dr ' yields 



where A+ := max{A,0}. Holder's inequality J /1/2 < ||/i||)9il|/2||p,7 applied to the functions /i 

' V'tv and f2 = <t>M 



(m— 1)t 



(55) 



Furthermore for any T g T/v and any z (^T, we obtain using ( 54 1 

1 



0^^(z)^^^.(z)^' = \T\^ ' ' {zf < ^||V(/ - 1^^ /)lli.(T)> 



hence 



TeTn 



E ^ /jiv(/-irV)iii.(T) = iiv(/-irrV)iri. 



(f2)- 



Elevating (55) to the power ^, injecting the above estimate and observing that = #7Ar < we 

thus obtain 

W^Nh^n) < l|V(/-l!;i;V)llL.(o)^"^- (56) 
Since the sequence {Tn)n>No is admissible, there exists a constant Ca > such that for all N and 
all T G T/v we have diam(r) < CaN^^. We introduce a subset of TJ(r C T/v which gathers the most 
degenerate triangles 

U^ITeTn; p{T) > uj{CaN---)^}, 



where oj is defined by (51 1. We denote by fl'j^ the portion of Q covered by T^- For all G il \ one 
has 

i^wiz) > L^,p{Trz) - Ci,\J uj{CaN^^). 
Hence, with the convention Lm,p{d"^ f{z)/ml) := Lm,p{'^z), 



II^ArllL-(n) 



> 



> 



\ to! 



C,\Iu(CaN-"2) 



L-(n\f2^) 



where C := max{L,„.p(7r2); z e 17} < 00. We next observe that — as +cx): indeed for all 

T € TIj we obtain using Lemma [lT] 



|T| < |req|diam(r)2p(r)-i < \T,^\CIN-'cu{CaN-2) 

i2 



cql 

. 1 X 1 



1 X 1 

+1 



Since #(7^) < N, we obtain < |roq|C^w(CAiV"5)^, which tends to as TV ^ cx). Therefore 



liminf ||-0/v|U-(n) > Jim 









\ to! y 











Combining this result with (56) we conclude the proof of the announced estimate (12). 
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3.2 A triangulation containing small periodic patches 

This subsection describes the construction of some triangulations by the aggregation of small periodic 
patches, which is a preliminary step for the proof of the estimate ( |13[ ) of Theorem[2]in the next subsection. 
The design of these triangulations is related to the anterior works IH [S] , yet it is adapted in order to keep 
under control the measure of silverness of the elements. The triangulations considered in this subsection 
are denoted by the letter V , instead of T, in order to avoid conflicts of notations with the next subsection. 

Our first lemma describes a family of meshes of a triangle R, which consist for the largest part of a 
periodic tiling based on a triangle T scaled down by a factor l/n, except for a few elements close to dR. 

Lemma 3.2. Let R and T he two triangles. There exists a family {VT,R.n)n>i, of conforming triangu- 
lations of R, and a constant Cr_t, such that the following holds. 

1. (Tights hounds on the cardinality ofVx.R.n md the diameter of its elements) 

#(Pt^^ ^ diani(r') < -(diam(T) +diam(i?)). (57) 

n-foo \T\ t'&Vt.r,^ n 

2. (Conformity) The vertices ofVT,R,n on the houndary of R are exactly those of the form ^a+(l— 
where < k < n and a, b are vertices of R. 

3. (Control of the houndary elements) Denote hy Rt.u C R the union of all the elements of VT,R,n 
which are not of the following form: the image of T hy a map of the form z t-^ Zq + az/n, for some 
zq € IR^, a e {-1, 1}. Then for all n > I 

\RT,n\<^^ and max S{T')<Cr^t- (58) 

Proof. See Appendix. □ 
We next introduce the concept of local shape specification, on the domain fl. 

Definition 3.3. A local shape specification is a (possibly discontinuous) map y Cz ^ f-^ Ty which asso- 
ciates a triangle Ty to each point y in the closure ofQ, and which satisfies the following properties. 

• The volume map ?/ G 51 > 171,1 G ]R*^ is continuous and positive. 

• The measure of degeneracy y G O h- > p{Ty) G -2?+ is uniformly bounded. 

The next proposition describes a sequence of triangulations adapted in a certain sense to a given local 
shape specification. This can be compared to the construction, evoked in §1.2[ of a triangulation adapted 



to a given Riemmannian metric in the sense of ( 35 ) . Given a point y G IR and compact set T we define 



df{(y,T) :— max{|y — z|; z G T}, which is the Haussdorf distance separating the sets {y} and T. 

Proposition 3.4. Let y Ty he a local shape specification. There exists a sequence {Vn)n>2 of trian- 
gulations of fl, a sequence {Sn)n>2 of positive reals converging to 0, and a constant Ca, satisfying the 
following properties. 

1. (Tight hounds on the cardinality ofVn and the diameter of its elements) 

lim ^^"^""^ = / ■— ^ and max diam(T) < — -. (59) 
n-^oo \Ty\ TeVr, n 

2. (Control of the houndary elements) Denote hy C 51 the union of all the elements T G Vn which 
are not of the following form: the image ofTy, for some ?/ G such that dfi{y,T) < Sn, by a map 
of the form z i— zq + az/n, for some zq G , a G {—1, 1}- Then for all n > 2 

|f2„| < — and max S(T) < Colnrt. 

n Tev„ 
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Proof. We consider a triangulation TZ^ of the polygonal domain f2 of minimal cardinality Nq. For each 
A: > 1 we denote by TZ'^ the triangulation of il of cardinality /c^iVo obtained by uniformly subdividing the 
elements of TZ^ into sub triangles. 

For each n > 1 we denot e by TZ'^ the triangulation of obtained as the union of the triangulations 



'PR,T,n described in Lemma 3.2 where R £ TZ and T = T^^ is the triangle specified by the local 



shape specification at the barycenter of R. Point 2 of Lemma 3.2 guarantees that this triangulation is 
conforming: there is no hanging node at the interfaces of the triangles R G TZ^ . Our next observation is 
that for each k,n> 1, using Lemma [lT] 



n max diam(r) < 3 majc (diam(i?) + diam(T^„)) < 3diam(f]) + sup {A/\T^q\)\Ty\p{Ty), 
which is finite and independent of k and n. For each k,n > 1 we define a real d{n, k) by the equality 



O \Ty\ 



(5(71, fc). 



Using (571 we obtain 



hm 6{n,k) = S{k) -.^ y M^ - [ ^ 



3.2 



for i? e 7^'= and T = T,^ , and 



Note that 5{k) — >■ as fc cx) since the map y ^ Ty is continuous. 

We denote by the union of the sets i?r,n described in Lemma 
by Cfe the sum of the corresponding constants Cr^t- We obtain using (158^ that for all n > 1 

\^n\ < — and max S{T) < Ck- 
n TeK^ 

We finally choose a sequence (fc(n))ra>2, such that fc(n) — >■ oo as n — >■ oo, and which increases "slowly" 
in the following sense: we require that for n sufficiently large one has 6{n, k{n)) < S{k{n)) + ^ and 
Cfc(n) !i Inn. Defining Vn '■= TZ^^^^ and (5„ := diam(f2)/fc(ri) we obtain the announced result. □ 



3.3 The upper error estimate (13) 

Throughout this section we consider a fixed real M > 1, and we introduce the collection of triangles 

Tm := {T- \T\ = 1, piT) < M, zt = 0} 

which is compact for the Haussdorf distance. We introduce a variant Lm of the shape function L^^p 
defined as follows: for each tt G Hm 

LmW min |1V(^ - l^-^ 7r)|U^(^). (60) 

-I t -l Af 

For any fixed tt e ]H,„ the map T i—^ ||V(7r — I™^^ t^)\\lp{t) is continuous with respect to the Haussdorff 
distance on the set of all triangles. Hence there exists a minimizing triangle, that we denote by T(7r) S 
Tj\/ , for the optimization problem appearing in ( 60 ) . 

Since all norms are equivalent on the finite dimensional space ]H„i, there exists a constant Cm such 
that for all tt e Hm 

sup ||V(7r - 7r)||iP(T) < Ca/||V^||. 

The function Lm is defined as the infimum of a family of C7\/-Lipschitz functions on H„i, hence is also 
CM-Lipschitz: \Lm{'^) — Lm{t^')\ < C'j\/||V7r — V7r'||. Furthermore for each 7r,7r' G H„,, one has since 
T(7r) e Tm 

l|V(^'-I?(;;7r')|U.(TW) < ||V(^-I™;;7r)|Up(^(,))+CM|k-7r'|| =LM(7r) + CM|l7r-7r'||. (61) 
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We consider the following local shape specification, see Definition |3.3[ 

y^Ty:= {LMiTTy) + M-^)-^^T{^y). 

In other words Ty is the isotropic scaling of the triangle TiiTy) by the factor [LMii^y) + M^^)^5. We 
thus have \Ty\ = (LMiTTy) + M~^)~'^ , which depends continuously on y G fl, and p{Ty) ~ p{T{ny)) < M. 

describes a sequence {Vn)n>2 of triangulations of VL attached to this local shape 



Proposition 



3.4 



specification, as well as a sequence ((5„)„>2 of positive numbers, a sequence {^n)n>2 of subdomains of £7, 
and a constant Ca- We recall that for all ti > 2 and all T E Vn 

diam(T) < Ca/n, S{T)<Calnn, |r2„| < Ca(lnn)/n. 

Using Point ii of Lemma |3.1| we obtain for each n > 1 and each triangle T e Vn 

l|V(/-irV)l|Lp(T) < aCr|T|?(lnn)/n™-i. 
Summing up the p-th power of the contributions of all the triangles T e Vn such that T C ri„, we obtain 



-,m— 1 ' 



and we observe that £„ — ;> as n oo. 

We now turn to the contribution of 51 \ f2„ to the error, and for that purpose we consider a triangle 
T which is the image of the triangle Ty, for some y € such that d-}i{y,T) < (5,1, by a map of the form 
z ^ zq + azjn, for some zq € H^, a G { — 1, !}• We have for any z £ T, using Point i of Lemma 



3.1 



liV(/ - I™-^ /)|Up(t) < ||V(^, - IJ^-' 7r,)|Up(T) + CMCa/n)\T\^M — 



where we used the invariance property (|8j) in the second line, and (61) in the last line. It follows that 
there exists an integer no such that ||V(/ — I™~^ /)||lj'(t) < for all n > no and all T EVn- Hence 

l|V(/ - r^;' /)|Up(o) < ||V(/ - I™-i /)|U.(a\n„) + ||V(/ - I^;' Dhnnr.) 



and therefore using ( 59 1 



limsup#(7'„) 



|V(/-i;rV)ll 



- er, n 



-(m-l) 



LP(0) 



< lim 



2 Sr, 



-,m— 1 



dy 

n \Ty\ 



(LMiTTy) + M-'Vdy 



(62) 



Observe that LMiTTy) converges decreasingly for each ?/ e f2 to Lm.pid™' Jiy)/m\) := L,„^p(7rj,) as 
M — >■ cxD. Given any fixed £ > 0, and using standard results on the convergence of integrals, we may 
therefore choose M = M(e) sufficiently large such that 



iLMi^y)+M-')^dy 



< 



d"^f 



■ e. 



(63) 



We denote by A'o the minimal cardinality of a triangulation of the polygonal domain 51, and we assume 
without loss of generality that #(7^2) = .^o- For each N > Nq, we denote by n(A^) the largest integer 
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such that #('Pn(Ar)) < N, and we set := Vn{N)- Observing that ^[T§)/N — > 1 as iV — >• oo, and 
combining (|62| with (63), we obtain the announced result (13): 



limsupTV— ||V(/-I^rV)llL.(o) =limsup#(P„(^^))— |1V(/-I™-^ f)\\Lr^n) 



< 



L 



+£. 



The admissibility of the sequence {T^)N>Na of triangulations immediately follows from (59). 



Conclusion 

In this paper, we have introduced asymptotic estimates for the finite element interpolation error measured 
in the W^'^ semi-norm, when the mesh is optimally adapted to a function of two variables and the degree 
of interpolation m — 1 is arbitrary. The approach used is an adaptation of the ideas developped in |18j 
for the LP interpolation error, and leads to asymptotically sharp error estimates, exposed in Theorems [T] 
and [2] These estimates involve a shape function Lm,p which generalises the determinant which appears 
in estimates for piecewise linear interpolation. The shape function has equivalents of polynomial form for 
all values of m, as established in §2.2| Up to a fixed multiplicative constant, our estimates can therefore 
be written under analytic form in terms of the derivatives of the function to be approximated. 

Future efforts will be devoted to the extension of these results to functions defined on a domain 
of dimension d > 2, which is partially done in Chapter 3 of [20] • Another challenge left open is the 
development of an anisotropic mesh generator with guarantees on the maximal angle of the elements, 
as evoked in |1.2| which could allow to apply the results of this paper in the context of adaptive mesh 
refinement for numerical simulations. 

APPENDIX 



A Proof of Lemma 13.2 



We construct the triangulation 7^T,i?,n of R as the union of three components: VT,R,n = S„ U £„ UT„. 
The elements of Bn cover the boundary of R, while the elements of I„ cover most of its interior and 
are included in i? \ Rt.u- The elements of play the role of a layer between I„ and Bn- Throughout 
this proof we denote by C = C{R,T) a generic constant independent of n, which may change from one 
occurrence to the next. 

We introduce the homothetic contraction i?„ of the triangle R by the factor 1 — n^^, with the same 
barycenter (i.e. the image of R by the map z z/? + (1 — n~^){z — zj{)). One easily checks that for all 
zedR 

LiamU ^"^(^'-^") ^ ^diam(j;), (64) 

where d{z, E) := inf{|2: - e|; e G i?} for any z £J\? and any i? C M^. 

Let zq be a vertex of the triangle T, and let u, v be two of its edge vectors. For each n > 1 we denote 
by In periodic tiling of built of the images of T by the maps z M- au + (3v + (zq + (t(z — zo))/n, 
where a, £ 7L and a S {—1, 1}. This tiling is built of translations of the triangle T, and of its symmetric 
with respect to the vertex zq, scaled by the factor 1/n. We define the collection of triangles 

In := {T' e II- T' C i?„}, 

which is illustrated in dark gray on Figure |3] Any T £ I„ satisfies diam(r') = diam(T)/n and S{T') = 
S(T). Furthermore these elements are of the form mentioned in Point 3 of the lemma, hence are included 
in i?\i?T,n- The elements of I„ cover the set {y £ R; n d{y, dR) > diam(T) + diam(i?)}, which area 



is larger than \R\ — C/n. We thus obtain the left part of (58). Observing that \T'\ = iTl/ri^ for each 
T' £ V, we obtain that #(I„)/n2 \R\/\T\ as n ^ oo. 
We define a collection Cn of convex polygons as follows: 

£° := {Rn fM' -r £ll and int(T') n dRn ^ 0}, 
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Figure 3: The triangles R and i?„ (left). The collections In (gray) of triangles and C'^ (white) of convex 
polygons (center left). The collection £„ of triangles is obtained by triangulating the elements of £^ 
(center right). The collection Bn of triangles covers R \ Rn (right). 

where int(i?) denotes the interior of a set C H^. The set £° is illustrated in white on Figure [s] (center 
left). The elements T' e Z" such that int(T') n dR^ ^ are included in the set G M^; nd{y,dR) < 
diam(i?) + diam(T)} which area is smaller than C/n, and their individual area is |r|/n^. Therefore 
#('C°) < Cn. The normals to faces of the elements of £° belong to a family of at most 6 elements: 
the normals to the faces of T, and to the faces of R. Therefore at most 6x5 different angles can 
appear in , and we denote the largest of these by a < TT. We denote by £„ the collection of triangles, 
illustrated Figure|3] (center right), obtained by triangulating each element of C'^^, which is a convex polygon 
with at most six faces. The angles of the triangles partitioning a convex polygon are smaller than the 
angles of this polygon, hence the angles of the elements of £„ are also bounded by a. Furthermore 
#('Cn) < 4#(£0) < Cn, and diam(r') < diam(r)/n for each T' G Cn- 

We denote by En the collection of n equidistributed points on each edge of R, described in Point 2 of 
Lemma [3. 2[ and we denote by E'^^ the collection of vertices of the triangles in £„ that fall on dR'^. For 
each point p e En, we draw an edge between p and the point of p' G E'^ which is the closest to p. Note 
that |p — p'l < d{p, Rn) + diam(T)/n < (diam(i?) + diam(r))/n. This produces a partition S° of \ i?„ 
into triangles and convex quadrilaterals, of diameter at most 

diam(i?)/n + 2(diam(i?) + diam(r))/n < 3(diam(i?) + diam(T))/n, 

since the distance between two consecutive points in En is at most diam(i?). We denote by Bn the 
collection of triangles, illustrated Figure [3] (right), obtained by triangulating each polygon K g Bn, of 
vertices Fn(£;„U£;;). We have #(S„) = #{En) + #iE'^) < 3n + 3#(/:^) < Cn. In order to conclude the 
proof of this lemma, we only need to show that the angles of the elements of Bn are uniformly bounded 
away from 77. 

We consider a triangle T' G Bn, we denote by L the length of the edge of T' included in dRUdRn, and 



by H the height of the triangle T' such that LH = 2\T'\. It follows from (64) that H > 2|T|/(3ndiam(T)) 



Let L' be another edge of T', and let 9 be the angle of T' between the edges L and L'. Then 

2\T'\ = LL' sine = LH, 

hence sin 6* > j,arnXT') — '^^^re c — 2|T|/(3n diam(T)) and C — 3(diam(i?) + diam(r)), which implies 
that arcsin(^) < < ir — arcsin(^). It follows that all the angles of T' are smaller than tt — arcsin(^) 
which concludes the proof. 



B Proof of Theorem 2.5 



We consider an arbitrary but fixed r > 2 and we define Sr ■— [|J + 1. It is established in Proposition 
2.1 of [IS] (equivalently Proposition 2.2.1 of 20J) that for any vr e the three following properties are 
equivalent 

" Krin) = 0, 

There exists a, /3 G IR and tt € H^-Sr such that tt = {ax + /??/)*'' tt, (65) 
There exists a sequence {4>n)n>o, 4>n G SL2, such that tt o 0„ — 0. 
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In addition the following invariance property is established in Theorem 2.6.3 of [20^: let Q be a polynomial 
on Mr such that Kr ^ \/\Q\ where d — deg Q. Then dr/2 is an integer and for all tt e Mr and all G M2 



Q(7ro0) = (det0)Tg(7r). 



(66) 



It follows that the polynomials {Qk)a<k<d, defined in (47), satisfy for all ^1,7:2 € Mr and all (f> G M2 



Qk{Tri o (f)^TT2 (j)) = (det(/)) 2 (5fc(7ri, 712). 
We define two functions on Mr x Mr 



if*(7ri,7r2) := 2d/ ^ Qfc(7ri,7r2)2 and i4:(7ri, vra) := '^Q(7rf + 7r|), 

V 0<fc<T 



(67) 



(68) 



where Q is an homogeneous polynomial on M2r such that K2r ^ \ Q, and d :— degQ- We show below 



that K ^ on Mr x Mr- Choosing r = to — 1 and combining this result with (44) concludes the proof 
of Theorer n 12.51 

Using (67) and remarking the invariance property Q{Tr o (j>) — (det (j))'^^Q{TT), for the same reasons as 



( 66 ) , we obtain for all 711,1^2 G Mr and all cf) e M2 



K{tti o 0, 7r2 o ( 
K.^ (tti o 0, 7r2 o ( 



= I det ( 
= det ( 



lK{'Ki,7T2), 
iK^{tTi,TT2) 



(69) 



Consider 7ri,7r2 G E^. The equality K{7Ti,tt2) = is equivalent to ttJ + Trf G M2r having a linear 
factor of multiplicity S2r = r + 1 (according to (65)), which is also equivalent to tti and 7:2 having a 
common linear factor of multiplicity Sr- 

On the other hand the equality K^{TTi,Tr2) = is equivalent to (5fc(7ri,7r2) — for all < fc < d. 



This is equivalent to Kr{ui:i + VK2) = for all u,w G IR (using (47)), which means that the polynomial 



MTTi + vt:2 G Mr has a linear factor of multiplicity s,- for all u,w G IR (using (65)). This is equivalent to 
TTi and 7^2 having a common linear factor of multiplicity s^. 
The following properties are therefore equivalent 

i^(^i,7r2) = 0, 

K^{tTi,TT2) = 0, 

There exists a, /? G R, and tti, 7r2 G Mr-s such that tti = [ax + PyY^'iri, and 7r2 = (ax + /?!/)'''■ 7r2. 

(70) 



Using (65) and the definition (68) of K, we find that these properties are also equivalent to 



if2,.(^? + 7r|) = 0, 

There exists a sequence (0„)ri>o, 4>n G SL2, such that (tti o 0„)^ + (7r2 o (/)„)^ — > 0, 
There exists a sequence {4>n)n>o, 4>n € SL2, such that tti o 0„ ^ and 7r2 o — >• 0. 



(71) 



We now define the norm ||(7ri,7r2)|| sup|,„|<]^ |(7ri(u), 7r2(w))| on H.^ x Mr, and the set 

J" := {(7ri,7r2) G x H,. ; || (tti, 7r2)|| = 1 and ||(7ri o 0,7r2 o (/))|| > 1 for aU <f) G SL2}. 
The set F is compact subset of Mr x H^, and K as well as if* do not vanish on F according to (70) and 



(71 ). Since these functions are continuous, there exists a constant Co > 1 such that 



Cn^K <K^< CqK on J". 



Let (7ri,7r2) G Mr x H^. If there exists a sequence 



(72) 
and 



t'n)n>0, 4>n & SL2, SUCh that TTi O 0„ 

7r2 o (/()„ — )■ 0, then X(7ri,7r2) = isr(0,0) = and K^{Tri,TT2) = K^{0,0) = using (69) and the continuity 
of K and if*. Otherwise, consider a sequence (0n)n>O7 € SL2, such that 



lim IKtti o 0„,7r2 o , 



inf IKtti o 0,7r2 o 0)||. 

pSfaL2 
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By compactness there exists a pair {ttijTt-z) G Hm x IH„i and a subsequence {(l>nk)k>o such that 

(tTi O (/)„^,7r2 O 0„J -> (7fl,7r2)- 



One easily checks that e J^. Using (|69|) we obtain 



= lini — 



K^.{'Ki,Tr2) n^oo ii:^(7ri o (/)„, 7r2 o 0„) i4r47ri,7r2) 



Using (72) and the homogeneity of K and X*, we obtain that Cq^K < < CqK on Mr x which 



concludes the proof. 
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